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Abstract: 

We analyse the kinematic regimes attainable for J/ 1\) and ip' production in p — A 
collisions. With this information, we specify the requirements for an experiment to 
study the interaction of physical charmonium states in nuclear matter, making use of a 
heavy ion beam incident on a hydrogen or deuterium target. 
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1. Introduction 

The production of charmonium states in hadr on- nucleus and nucleus- nucleus colli- 
sions provides an excellent tool to study different aspects of strongly interacting matter. 
At high energy, charmonium production proceeds through gluon fusion to a colour octet 
cc state (Fig. la), which subsequently neutralises its colour and obtains the right quan- 
tum numbers by interacting with the colour field of the collision [1] . At lower energies, 
quark-antiquark annihilation adds significant contributions (Fig. lb), which for large 
\xf\ become dominant. In both parton processes, the first production stage is a colour 
octet cc state which needs a certain time before it becomes colour neutral. Consider 
p — A collisions leading to the production of charmonia fast in the rest frame of the target 
nucleus. The nuclear medium then sees only the passage of a fast coloured cc pair, with 
which it interacts as with any colour charge. There is nothing charmonium-specific 
about this interaction, and in particular, the nucleus does not know what hadronic 
state this pair will turn into later on. To study the effect of the medium on the dif- 
ferent physical charmonium states, these have to be sufficiently slow to be fully formed 
inside the medium. In the first part of this note, we will identify the kinematic regions 
corresponding to these different situations. As we shall show, all hadron-nucleus colli- 
sions studied so far provide information only about the passage of the coloured cc pair 
through nuclear matter; nothing is known experimentally about the interaction of fully 
formed charmonia with the nuclear medium. 

This lack of information is particularly serious if one wants to use charmonium 
production as a probe for deconfinement in nucleus-nucleus collisions. It was predicted 
[2] that quark-gluon plasma formation leads to J/ip suppression in nuclear collisions, 
and such a suppression was subsequently in fact observed [3]. This observation in turn 
triggered a "conventional" explanation [4], in which the suppression was attributed 
to absorption of the J/ip in a dense hadronic (and hence confined) medium. More 
recently it was shown [5] that the theory of heavy quarkonium interactions with light 
hadrons leads to a break-up cross section which rules out such an absorption. It was also 
noted, however, that a direct experimental test of the crucial cross section behaviour 
is still lacking. The theory on which it is based becomes exact in the heavy quark 
limit, and although some related checks were favorable, it is certainly necessary to 
verify experimentally that the charm quark mass is already sufficiently heavy. In the 
second part of this note, we compare the predictions of heavy quark QCD with previous 
geometric estimates, in which the QCD threshold behaviour had not been taken into 
account. From this comparison it becomes clear that an experimental verification of the 
predicted J/^-nucleon cross section behaviour is possible. After a discussion of possible 
initial state modifications (EMC effect) we then specify the kinematic range and further 
aspects of an experiment specifically devoted to this question. 
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2. Production Kinematics 



The time needed by the colour octet cc state to become colour neutral can be 
estimated in terms of the energy of a further gluon "evaporated" to produce a colour 
singlet state [6] . In the rest frame of the cc, this time r becomes 

to = -^=; (l) 

m c is the charm quark mass and u) the gluon energy. In confined matter, u> > Aqcd, 
making the colour neutralisation time in the rest frame of the cc pair tq — 0.25 fm. In 
the rest frame of the nucleus, the cc travels in the time To the distance 

do = (§) t , (2) 

with Pa denoting the momentum and M the mass of the cc state. From Eq. (2) it is thus 
clear that sufficiently fast cc are still coloured when they leave the nuclear environment. 
The Feynman variable Si? of the charmonium state and hence of the corresponding 

cc is 

xf = , (3) 



where 

P = lP A - l(3yJP 2 A + M* (4) 

is the center of mass (cms) momentum of the cc; 7 and (3 specify the transformation 
from the nuclear rest frame to the cms. The maximum cms momentum the cc can have 
is determined by the relation 



* = (y/PLx + M 2 + ^Pl ax + 4m 2 ) 2 , (5) 

with s for the squared incident cms energy and m for the nucleon mass; Eq. (5) corre- 
sponds to the configuration of both nucleons going forward and the cc backward. 

We shall now consider the specific case of a 160 GeV/c proton beam incident on a 
nuclear target, corresponding to = 17.4 GeV. At this energy, it is possible to study 
p — A collisions at the CERN-SPS with a nuclear beam incident on a hydrogen target 
as well as the conventional inverse. In Table 1, we list the resulting windows for the 
production of J/ip, Xc and ip'. We thus obtain 8.3 GeV as the maximum momentum 
which a J ftp can have in the cms; this means that in the nuclear rest frame, it has a 
minimum momentum of 4.7 GeV. In terms of rapidity, the J/ip production range at this 
energy is |y C ms| < 1-7; in the region between |y C ms| = 1-7 and the maximum rapidity 
|2/cms| = 2.9, J/tjj production is kinematically not possible. 

Since we want to study the interaction of physical charmonium states with the 
nuclear medium, we want to avoid the kinematic regime in which the cc can interact 
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with the nucleus already in its colour octet phase. Hence we require d < 1.5 fm, so 
that the cc has become colourless before it leaves the range of the nucleon on which 
it was produced. Using the above relations, this implies at the noted collision energy 
Xi? < - 0.2; for all larger xf, at least a partial colour interaction with the nuclear 
medium remains. In the region —0.2 < xf < 0, this is effectively negligible, since the 
distance the cc travels as coloured state is only about 2 fm at if = 0. The colour 
phase increases rapidly with increasing xj?, and for ip> 0.2, the cc is coloured on its 
entire path through the nucleus [6]. Hence the nuclear effects on J/ip and ip' observed 
in this region must be identical. As already mentioned, all presently available data [7,8] 
on charmonium production in hadron-nucleus collisons fall into this region, and they in 
fact observe the same nuclear suppression of J/ip and ip' production. 

If we wish to study the interaction of fully formed physical resonances with a 
nuclear medium, it is clear then that we must require if < —0.2; this is a necessary, 
not a sufficient condition. For the latter, we have to assure that the cc has not only 
lost its colour, but that it has also had sufficient time to attain its full physical size. 
An estimate of the resonance formation time r r for the J/ifr is obtained from potential 
theory [9], with T r (J/ip) ~ 0.35 fm; this agrees with the rather model-independent upper 
bound based on the uncertainty relation and the separation between ground state and 
first excited state. The Xc and ip' values are larger; the noted potential theory studies 
are in accord with about 1 fm for both. With these values and the relation 



we find that the J/ ip passes through the entire nuclear medium as a fully formed reso- 
nance for xf< — 0.45, but that at the energy under consideration the Xc and tp' never 
reach that stage. Although they become colourless in the medium when xp < —0.2, they 
have not quite reached full physical size when leaving the medium even for Xf ~ —1. 

3. Charmonium Absorption 

The effect of the nuclear medium on the passing coloured cc state has been measured 
at two different energies [7,8]; it can be understood in terms of the energy loss of the 
colour octet in a confined medium [6] together with quantum-mechanical coherence 
effects (nuclear shadowing) [10]. The results are schematically illustrated in Fig. 2 in 
the range xf > 0.1; as noted, the effects on J/ip , Xc and ip' are identical. Below 
xf = 0, the cc passes the nucleus in part as a colourless state, and below xf — —0.2, 
it is colourless on its entire path through the nucleus. We can parametrise the survival 
probability for state i (i = J/ip, Xo as 



where no = 0.17 fm 3 is the normal nuclear density, L = (3/4)Ra = (3/4)1.15 A 1 / 3 
the average path length in the nucleus and Oi the absorption cross section of the state 




(6) 



Si = exp{-n <JiL}, 



(7) 
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% in nuclear matter.* If the charmonium state is fully formed before it leaves the range 
of the nucleon at which it was produced, <7j is simply the charmonium-nucleon cross 
section. If it is not yet fully formed, the effective cross section will be smaller, vanishing 
in the colour transparency limit of a pointlike colour singlet cc state. In principle, the 
interaction of this evolving resonance should be studied quantum-mechanically. Here 
we shall for simplicity parametrise the cross section as function of the distance d which 
the state has travelled [12,13], so that 

where Oi is the fully developed cross section and L = [(Pa/M^Tt — 1] the effective 
distance travelled until full resonance formation. Eq. (8) holds for d < L; for d > L, 
(Ji(d) = Oi. Using this parametrisation, Eq. (7) is replaced by 

2 _ 

Si = exp{-n a t [L - -L t ]} (9) 

whenever d < L. To determine the actual survival probabilities, we now still need the 
cross sections Uj for the fully formed resonances colliding with nucleons. 

In the geometric approach to charmonium absorption, is assumed to be the 
total high energy collision cross section. This can be estimated by geometric arguments 
[17], giving ctj/tp ~ 2.5 mb, a Xc ~ 6.0 mb and a^> ~ 9.2 mb. With this, the picture is 
complete: for each value of xf, we have the corresponding momentum Pa in the nuclear 
rest frame; from that we get in turn the distance travelled in the medium as nascent 
resonance and the associated cross section. The difference in asymptotic cross sections 
and the different formation times then lead to the absorption patterns shown in Fig. 2 
for negative xp. The region around xp = is drawn as continous interpolation of the 
forms at positive and negative xp. 

The crucial assumption leading to the behaviour just discussed is that the cross sec- 
tions for the interaction of charmonium states with nucleons attain their full asymptotic 
values at threshold. The invariant energy ^/I for the quarkonium-nucleon interaction is 
below 8-9 GeV for the region of xp < 0; for xp < —0.5, < 6 GeV. The threshold 
for DD production in charmonium-nucleon interactions is 4.7 GeV. Hence it is the in- 
elastic cross section near threshold that matters for the survival probability at negative 
xp. For the J/ip, and perhaps also for the Xc, this cross section is calculable by short 
distance QCD [14-16] [5]. A break-up requires the interaction with hard gluons, but 

* We note here that if the form (7) is used to parametrise the nuclear suppression of 
charmonium production at positive xp [11], then Oi becomes the i-independent cross 
section for the interaction of a colour octet cc with a nucleon. This can be quite large, 
but it is not related to that for the interaction of physical quarkonium states with 
nucleons. 
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gluons confined to slow hadrons in the charmonium rest frame are generally very soft. 
As a consequence, the break-up cross sections near threshold are expected to be much 
smaller than their asymptotic values. For the J/ip, the heavy quark theory predicts [5] 



oj/4,n{s) ~ 2.5 mb x 1 - 



. 6.5 



(10) 



where ejy^ = 2Me> — Mj/^ ~ 0.64 GeV is the J /if) binding energy and m the nucleon 
mass. The behaviour of this cross section is shown in Fig. 3 as function of the momentum 
Pjv of a nucleon incident on a J/ip at rest. For — 4 GeV/c, corresponding to y 7 ! = 6 
GeV, it is almost two orders of magnitude below its asymptotic value. For a cross section 
of this size, the survival probability is in good approximation unity, and that is the basis 
of the prediction that confined matter is transparent for J ftp's. 

The theoretical basis for the use of short distance QCD is much less reliable for 
the Xci since here the binding energy is just around Aqcd- Keeping this in mind, we 
shall nevertheless see what a corresponding analysis leads to. Instead of Eq. (10) we 
now obtain 



<j x n(s) ~ 11.3 mb x ( 1 — 



2M x (m + e x m 6 - 5 

{s-MD \) ' W 

where e x = 2M D — M x ~ 0.24 GeV now denotes the binding energy of the Xc- The 
asymptotic value is a factor two larger than the geometric estimate; this is a consequence 
of the fact that short distance QCD [14,15] leads to higher powers in the bound state 
radii than just r 2 . The behaviour of the XcN cross section (11) is also shown in Fig. 3. 

The ift' state lies essentially at the open charm threshold and hence has a binding 
energy much less than Aqcd! it is therefore definitely not calculable in short distance 
QCD. The asymptotic form ifj'^ shown in Fig. 2 may thus be a reasonable estimate 
here. 

The resulting survival probabilies Sj/^xp) and S x {xp) are compared in Fig. 2 to 
those obtained above in the geometric approach. In the region around xp = we again 
make an estimate taking into account both the vanishing of the colour interaction with 
decreasing xf and the growing charmonium- nucleon cross section with increasing xf- 

4. EMC Effect 

The survival probabilities shown in Fig. 2 are related to the measured p — A and 
p — p production cross sections through 

(daf A /dx F ) _ g A (x 2 ) , (r\q(,r\ C\9\ 

A(daf/dx F ) ~ ~^) S%{XF) = ^M* 2 )^)' ( 12 ) 

where (7a(^2) and g p (x2) are the parton distribution functions in nuclear and proton 
target, respectively. The fractional parton momentum X2 is given by 



*2 = -(^(xl + (4Mi/s))±x F ) (13) 
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in terms of the variables xp and s; the plus sign holds for positive, the minus sign 
for negative xp. At the energy under discussion above (y/s = 17.4 GeV) and for 
%f < 0, we have xi > 0.15; we are thus above the region in which quantum-mechanical 
coherence effects (nuclear shadowing or antishadowing) play a role [10]. However, we 
know from deep inelastic scattering on nuclear targets that for x<i > 0.15 the quark 
parton distributions in nuclei are modified in comparison to those in nucleons. This 
modification, generally denoted as EMC effect [18], is shown in Fig. 4. 

Although the EMC effect has so far been observed only for quarks, it is to be 
expected that gluons will exhibit a similar behaviour, so that the initial state factor 
Ra/ p {. x 2) will introduce a further xp variation in addition to that coming from Si{xp). 
With increasing \xf\, charmonium production is more and more due to quark-antiquark 
annihilation rather than to gluon fusion; the two contributions become approximately 
equal around \xp\ = 0.5, and as \xp\ — > 1, the qq contribution is dominant [1,19]. We 
shall here assume that quark and gluon distributions behave in the same way and use 
the quark form of Ra/ p m the whole region —1.0 < xp < 0; this gives us a prediction 
for the behaviour of the cross section ratio in that region. It is illustrated in Fig. 5 for 
the directly produced J ftp state in the two cases considered. 

As a cross check of the form of Ra/ p (x2), the EMC effect in this region can also 
be studied independently. Measuring Drell-Yan dilepton production there provides 
Ra/p( x 2) for quarks directly, without any final state modification. A measurement 
of open charm production leads to Ra/ p ( x 2) m just the same superposition of quark- 
antiquark annihilation and gluon fusion as in charmonium production, but again without 
any final state effect. Separate measurements of Drell-Yan and/or open charm produc- 
tion would thus determine the EMC modification without addititional final state effects. 
Such measurements can therefore be used to remove the EMC modification of J/ip and 
ifj' production data, which can then be compared directly to the predictions shown in 
Fig. 2. 

Such measurements would moreover be of considerable interest in themselves. At 
positive xf, little or no nuclear modifications are seen in Drell-Yan production [20]. The 
same experiment also showed no nuclear effects on open charm production. These two 
facts are part of the empirical basis for the claim [10] that the considerable suppression of 
charmonium production in p — A collisions at positive xp is due to quantum- mechanical 
coherence effects and energy loss suffered by the passing virtual colour charge, and not 
to factorisable parton distribution function changes. The measurement of factorisable 
modifications for negative x p would provide further support to this interpretation. 

5. Experimental Aspects 

The study of charmonia or Drell-Yan dileptons slow in the nuclear rest frame has 
up to now been essentially impossible. Both require the detection of slow dileptons, 
and for this the abundance of slow hadrons constitutes an overwhelming background. 
In the case of fast dileptons, a hadron absorber can eliminate these, and hence all p — A 
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studies were restricted to dilepton pairs of more than 20 GeV in the rest frame of the 
nuclear target. This in turn restricts the available data to xp> 0. 

The advent of the P6-beam at the CERN-SPS has removed this constraint. With 
the P6-beam incident on a hydrogen (or deuterium) target, the nuclear rest frame moves 
with a lab momentum of 160 GeV. Hence now charmonia and their decay dileptons are 
very fast in the lab system and will thus pass the hadron absorber. The window for 
such measurements is evident from Table 1, and the expected behaviour for positive 
xp is that shown in Figs. 2 for negative xp. For clarity, we denote the variables in the 
case of a A-beam incident on a hydrogen target by a superscript A. Thus the region 
of greatest interest is 0.4 > x F > 1; in terms of rapidity, this corresponds to the range 
3.6 > jfa > 4.5. 

The above predictions for J/ifj production correspond to directly produced IS cc 
resonances. The J/tj; peak observed in the measured dilepton spectrum is about 60% 
due to this origin and about 40% due to direct Xc production with subsequent decay 
X c — > J/ip + 7 (see [1]). Unless it is possible to measure Xc production independently 
[21], the data will contain a superposition of direct 15 J/ip production and J/ip 's from 
Xc decay.* To estimate the effect of this in the two scenarios considered here, we simply 
add the corresponding predictions with the noted 60/40 weights; the result is shown in 
Fig. 6. As seen, for xp < —0.4, the two approaches differ qualitatively in their functional 
form and quantitatively by more than 20%. An experimental test should therefore be 
possible. 

Finally we note that the analysis proposed here is a comparison of production data 
from a heavy nuclear beam on a hydrogen or deuterium target with that from a proton 
beam on the same target. Although there are measurements for the latter (see [1] for 
a compilation), it seems very desirable to obtain both A and p beam data in the same 
experiment, in order to avoid acceptance uncertainties. 

6. Conclusions 

The study of charmonium production at low momenta in the nuclear rest frame, 
together with that of Drell-Yan dileptons and open charm, opens a completely unex- 
plored region of the behaviour of hard probes in nuclear matter. Such studies have 
become experimentaly feasible only with the advent of the P6-beam at the CERN-SPS, 
and since they require very forward measurements, they would not be easy to carry out 
in future collider experiments. The results of such a program would have a decisive 
impact on at least three different topics: 

- They would test directly the heavy quark theory prediction of strong threshold 

damping for the interaction of J/ip : s with light hadrons. This is crucial for the use 

of charmonium production as deconfinement probe. 



* In addition, a few percent (~ 5 %) will come from ip' decays. We neglect this in 
both scenarios. 
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- They would allow a study of the EMC effect by cc production and thus provide a 
complementary tool to deep inelastic scattering on nuclear targets. This is impor- 
tant for the investigation of nuclear modifications of parton distributions. 

- By comparing nuclear effects on charmonia with those on Drell-Yan production, 
the role of modified parton distributions could be tested directly. Fast cc pairs (in 
the nuclear rest frame) show a very different behaviour than fast dileptons, and 
this rules out the factorisable modification of parton distributions there. For slow 
cc and dileptons, the behaviour is predicted to be similar and and the modifications 
factorisable. 

The experimental investigation of Pb — p collisions could thus do much to further our 
understanding of the effect of a confined (nuclear) medium on the parton structure of 
hadronic systems. 
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Table 1: 

Windows for Charmonium Production in p — A Collisions at ^fs = 17.4 GeV 



Variable 




Xc 




rail [GeV] 


8.31 


8.24 


8.19 


|,,max 1 

lycMsl 


1.71 


1.59 


1.54 


Pf n [GeV] 


4.70 


6.12 


6.89 


„ min 

Va 


1.21 


1.33 


1.38 


.max 
S/lab 


4.63 


4.51 


4.45 
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Figure Captions 

Fig. 1: J/ip production by gluon fusion (a) and by quark-antiquark annihilation (b). 

Fig. 2: Charmonium production in p — Pb collisions at 160 GeV incident beam energy. 
The suppression for xf > is the same for J/i/j, if)' and Xc] data are for J/ip production 
at 200 GeV beam energy [7]. The predicted suppression for xp < is shown for 
asymptotic cross sections (solid lines) and for cross sections from short distance QCD 
(dashed lines). 

Fig. 3: Dissociation cross sections for J/^-nucleon and x c -nucleon interactions as func- 
tion of the momentum Pn of a nucleon incident on a charmonium at rest. 

Fig. 4: EMC effect: ratio Ra/ p = 9A(x)/g p (x) of quark distribution function gA(x) in 
heavy nuclei to that in nucleons g p (x) [18]. 

Fig. 5: Direct J/ip production in p — Pb collisions in the geometric approach and from 
short distance QCD, with (solid line) and without (dashed line) EMC effect modifica- 
tion. 

Fig. 6: J /if) production in p — Pb collisions in the geometric scenario and in short 
distance QCD, with 40% of the observed J/ip coming from Xc decay; without EMC 
effect. 
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